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1. Summary 

Cell damage by high linear energy transfer (LET) 
radiations has been described by a phenomenolog- 
ical model (track theory) for more than 20 years. 
With track theory, molecules of biological signif- 
icance (dry enzymes and viruses) act as one-hit 
detectors. Recent additions to the class of one- 
hit detectors are Escherichia coli B , and single- 
and double-strand breaks in SV-40 virus in EO 
buffer, where indirect effects predominate. The re- 
sponse of cells (survival, transformation, and chro- 
mosome aberration) to these radiations is typically 
described by a four-parameter model whose numer- 
ical values are determined from the equations of 
the theory being fit to experimental data at high 
dose (typically above 1 Gy), with the cells bom- 
barded by gamma rays and high-charge-and-cnergy 
(HZE) particle beams, of the widest possible dynamic 
range. Once these parameters are determined, the 
model predicts cellular response in any radiation en- 
vironment for which the particle-energy spectrum is 
known. The important feature of this track structure 
model is its ability to estimate from a limited set of 
laboratory data the response of a complex radiation 
environment with many components. For example, 
we have calculated cell survival after neutron irra- 
diation with mixtures of neutrons and gamma rays 
and cell survival and transformation after irradiation 
with HZE of different energies. The model does not, 
yet include cellular repair, although some hopeful ap- 
proaches to repair dependence are now being devel- 
oped. It also docs not include cancer induction be- 
cause the available data give neither the number of 
cells at risk nor the number of cancers induced and 
are thus not suited to our formulation. 

Most recently, NASA Langley models of HZE 
beams, including projectile and target fragmentation 
types, have been joined with the track structure 
model. This combination has been tested with good 
success against ground-based radiobiological data for 
cell survival after irradiation with protons and HZE 
beams. And whereas our earlier beam model failed 
downstream of the Bragg peak (for both protons 
and heavy ions) for want of a proper description of 
fragmentation, the current model succeeds. 

Based on this experimental validation of our pro- 
cedures, we have initiated calculations of cellular 
damage in spaceflight from solar protons and galac- 
tic cosmic rays. In this paper we incorporate NASA 
Langley models of cosmic rays, beam penetration, 
and projectile and target fragmentation with track 
theory. The essential radiobiological theme is that 
knowledge of parameters found at high doses en- 
ables us to calculate the response of cells at the low- 


est possible doses of HZE particles when only intra- 
track (ion-kill) effects arc involved. Our procedures 
here too have ground-based experimental validation, 
wherein measurements made with protons and alphas 
of relative biological effectiveness (RBE) of the sur- 
vival of C3H10TV2 cells, at doses down to 0.01 Gy, 
were consistent with our predictions that were based 
on survival measurements made at high doses with 
gamma rays and HZE’s. 

2. Introduction 

Detectors of radiation differ according to whether 
single-particle response is normally observed, as with 
nuclear emulsions, solid-state nuclear track detectors, 
and scintillation counters, or whether the response is 
to beams of particles or photons in a gross macro- 
scopic irradiation, as in radiobiology or in the alter- 
ation of bulk material properties by radiation. In 
the former case it is more natural to think in terms 
of track structure, whereas in the latter case one fre- 
quently refers to macroscopic dose (Katz 1978b). Re- 
sponse is then correlated to the physical description 
of these stimuli. It is common to relate response 
to energy deposition (dose). Problems arise because 
response depends not only on total energy deposi- 
tion but also on the microscopic structure of that 
deposition and its time development. One analysis 
of these details is called microdosimetry , the study 
of energy deposition in small volumes, which stim- 
ulated many investigations. An alternate procedure 
favored in this paper relates the observed effect to 
track structure for individual particles, and the effect 
then may be related to macroscopic dose for gross 
irradiations. These perspectives are principally re- 
ported in the several Symposia on Microdosimetry 
sponsored by the Commission of European Commu- 
nities (Katz and Huang 1991; Katz, Dunn, and Sin- 
clair 1985). 

The galactic cosmic ray (GCR) environment is 
the most complicated mixture of radiation compo- 
nents known. It is doubtful that the GCR environ- 
ment will ever be adequately simulated in the labo- 
ratory for biological experiments. The primary role 
of track-structure models will be to extrapolate lab- 
oratory response data to the GCR environment for 
the estimation of risk to biological tissues in space 
exposure. We believe this extrapolation will be a 
more practical approach to the problem of additivity 
of response of disparate components than the usual 
quality- fact or approach based on relative biological 
effectiveness (RBE), which has been used with lim- 
ited success in terrestrial radiation protection. 

This paper is organized as follows: We first 

discuss several of the approaches for describing 



radiation insult and emphasize the question of addi- 
tivity of radiation components. The track-structure 
model and its mathematical formalism are then de- 
scribed. The Langley models of high-charge- and- 
energy (HZE) transport are compared with cell sur- 
vival experiments and predictions of cell damage in 
space flight from galactic cosmic rays are discussed. 
Finally, significant conclusions are given. 

Research at the University of Nebraska is sup- 
ported by the U.S. Department of Energy. 

3. Symbols and Abbreviations 

A mass number 

«o sensitive site radius, cm 

C hit ted ness 

D absorbed dose, Gy 

D average radial dose, Gy 

D r X-ray dose, Gy 

D~ gam in a- ray dose, Gy 

E energy, MeV 

E { } radiosensitivity parameter, Gy 

E[). average energy of fragment j, MeV 

F fluence, cm 2 

/ secondary ion spectrum, MeV -1 

G number of observed events per unit of 

energy deposited 

GCR galactic cosmic ray 

HZE high charge and energy 

L stopping power, MeV/cm 

LET linear energy transfer, MeV/cm 

m target number 

N cell population after exposure 

Nq initial cell population 

n number of observable events per 

centimeter of path length 

P ion-kill probability 

p proton 

QF quality factor 

RBE relative biological effectiveness 

RBE ext extrapolated relative biological 

effectiveness 


S target cross section, cm 2 

So plateau value of the target cross 

section 

TLD thermoluminance detector 

t radial distance, cm 

x position, cm 

Z charge number 

Z* effective charge number 

j3 relative particle velocity 

k track structure parameter 

n z ion-kill fraction 

n 7 gamma-kill fraction 

E macroscopic nuclear interaction cross 

section, cm -1 

a action cross section, cm 2 

o 

(7 e xt extrapolated action cross section, cm 

<tq plateau value of action cross section, 

cm 2 

a* effective action cross section, cm 2 

r mission duration, yr 

particle flux spectrum, #/MeV-cm 2 -s 

Subscripts: 

j fragment label 

n neutron 

4. Concepts From Radiobiology 

4.1 Microdosimetry: Energy Deposition 
in Small Volumes 

One way to analyze stimuli to biological systems 
is to examine the details of energy deposition in small 
volumes that are sized to represent what are thought 
to be critical targets within the cell. Experimentally 
small, gaseous proportional counters are used whose 
diameters, scaled to the density of tissue, range from 
micrometers to nanometers in unit-density material. 
The critical targets are then considered to be the 
nucleus of a mammalian cell, a chromosome, or a 
small region of DNA. The fluctuations of the energy 
deposited within the small target region are assumed 
to be related to biological response. A Monte Carlo 
simulation of a radiation field can yield a similar 
decomposition. 
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Even when one has a complete microdosimetric 
description of the radiation environment, the prob- 
lem remains as to how that description may be inter- 
preted to predict the response of a detector. As yet 
we have no means of calibrating response in terms 
of the statistical distribution of energy depositions 
in small volumes; nor do we know what volume is 
appropriate. It is on this level that microdosimetry 
has not been able to make extensive quantitative pre- 
dictions, nor has it been able to yield calculations of 
cross section. But microdosimetry has yielded many 
interesting insights into the structure of a radiation 
field. (See, for example, Goodhead 1988.) The small 
counter has found applications in monitoring neutron 
beams that are used in radiotherapy and in other ra- 
diation fields, including spacecraft and high-altitude 
aircraft. Most instruments are used in practice 
to derive averages over quality factors discussed in 
section 2.5. 

4.2 Cross Section 

A second approach to analyze stimuli to biologi- 
cal systems is to attempt to mimic the kind of logical 
structure that is used in experiments in physics, that 
is, to describe the relevant interactions through the 
concept of an interactive cross section. We imag- 
ine that a projectile passing down a channel 1 cm 2 in 
area interacts with a target located somewhere within 
that channel, and we then measure the fraction of 
successes after a large number of identical repeated 
trials. The probability of success is represented as 
though it is a geometrical target, as the cross sec- 
tion a in square centimeters to the cross-sectional 
area of the channel. We then speak of the action cross 
section even if the observed end point is achieved as 
a result of many internal changes stimulated by the 
initial interaction. We make no attempt to exam- 
ine the internal processes mechanistically: the target 
is a black box. We know only the incident radia- 
tion and the observed end point. In radiobiology the 
concept of action cross section is sometimes used in 
ways that depart from its original physical meaning. 
This difference can lead to misinterpretations of ex- 
perimental data (Katz 1990). Curtis et al. (1992) 
have recommended an additivity formalism based on 
a limited set of data for harderian gland tumorige- 
nesis by using a cross-section-like formalism as an 
alternative to the use of RBE. 

4.3 G - Value 

When a projectile impinges on a thin slice of 
matter containing JVq targets per cubic centimeter, 
the number of observable events per centimeter of 
path length is n = aN . (If the stopping power is 


L, the number of observed events per unit of energy 
deposited, the G-value, is G = n/L — aN/L. The 
cross section is a function of the medium, the end 
point, and the character of the projectile as follows: 
if a photon its energy, if a naked charged particle its 
charge and speed, and if a nucleus partially clothed 
with electrons its effective charge and speed. This 
formulation of the G-value lias been used in the 
analysis of heavy ion radiolvsis (Katz and Huang 
1989). In dealing with liquids, where the meaning 
of N may be obscure, we have calculated the G- 
value for heavy ion bombardment from calculated 
values of the RBE and known G- values for gamma 
irradiation, as in the Fricke dosimeter (Katz, Sinclair, 
and Waligorski 1986). In other cases we have' tried 
to relate N to a fitted target size. Nuclear collisions 
are here neglected except as a source of charged 
fragments. 

If the G-value is normalized to molecular weight 
and expressed as events per radian per dalton rather 
than as events per 100 eV, we find it proportional 
to the RBE for dry one-hit detectors for which the 
target molecular weight is equal to the true molecular 
weight (Katz 1991). 

4.4 Relative Biological Effectiveness and 

Radiation Quality Factor 

When intense neutron environments became 
available for study, the existing body of biological 
response data was mainly for X-ray and gamma-ray 
sources. The first efforts at protection attempted to 
scale the known X-ray and gamma-ray risks accord- 
ing to equivalent neutron dose, and this scaling gave 
rise to the concepts of RBE and radiation quality. 
For the case of space radiations, we should like to re- 
late the response of our detectors to energetic heavy 
ions (high LET radiations) to their response to pho- 
tons and electrons (low LET radiations). In radio- 
biology the ratio of the dose of gamma rays to that of 
another radiation that produced the same observed 
end point is called the relative biological effectiveness 
(RBE). In radiotherapy this quantity is frequently 
taken to be a property of the two radiation fields, but 
it depends on the dose level, the dose rate, and the 
end point as well. An extension of this idea used in 
radiation protection is called the quality factor (QF), 
which is an estimate of the upper limit of RBE val- 
ues for a selected set of biological end points judged 
relevant to human risk and taken solely as a func- 
tion of LET. An important, unresolved question is 
whether an upper limit or maximum RBE is achieved 
at the low exposures of interest, for radiation protec- 
tion. The QF is used to convert a measured dose in 
gray into an effective dose reported (not measured) 
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in sieverts. A quantity in sieverts is not directly mea- 
surable and thus violates the underlying philosophy 
through which physical units are defined. The re- 
deeming quality of the sievert is that risk estimates 
based on the unit should be conservative, providing 
QF is adequately defined. The conservative nature 
of the method may also be an unacceptable burden 
in many operations, especially in space where it may 
result in severe penalties in the form of substantial 
increases in shielding. 

4.5 Radiation Dose and Quality 

An irradiation with photons leads to secondary 
electrons randomly dispersed through a medium. 
The initial energy spectrum of these electrons and 
their path lengths depend on the initial energy spec- 
trum of the photons. An irradiation with a beam 
of heavy ions yields a random distribution of heavy 
ion paths, with the secondary electrons (delta rays) 
clustered around each ion's path (correlation effects) 
and having a different energy distribution. Hence 
the delta rays are not truly randomly distributed, 
but are clustered about the paths of heavy ions. 
The difference in the spectrum of secondary elec- 
trons from photons and ions leads to a basic dif- 
ference in the manner in which their effects are ap- 
proached statistically. In the track structure model, 
gamma kill describes the effects of a random distri- 
bution of spatially uncorrelated electrons, and ion 
kill describes the effects of spatial correlations within 
single-particle tracks. At high fluences of low LET 
ions, where only a fraction of the intersected targets 
is inactivated, we attribute the effect of the sparsely 
distributed and overlapping delta rays from several 
ions to randomly distributed electrons and consider 
gamma kill to be responsible for part of the effect 
from beams of some heavy ions. 

During irradiation with photons, minutes may 
elapse before secondary electrons traverse from dif- 
ferent photons through a target (uncorrelated tem- 
poral events). During irradiation with heavy ions, 
however, a single ion and its delta rays pass through 
a target in an extremely short time because the pro- 
jectile moves at nearly the speed of light through a 
target whose diameter is on the order of 1 fim. These 
differences in time separation lead to a variation in 
target response to radiations of different admixtures 
of photons and heavy ions at the same dose. When 
the temporal correlation time is on the order of the 
cell repair time, then response also varies with dose 
rate or fractionation schedule. The dependence on 
dose rate or fractionation schedule is due in part to 
radiation quality. The QF is taken as the low-dose 
rate limit of RBE’s in an attempt to normalize the 


biological effects of radiations of different qualities. 
Choosing the quality factor in this manner thus as- 
sumes that it is logical to represent the response of 
a detector as a product of two separate factors, dose 
and quality. 

4.6 Radial Dose Distribution 

To make track-structure calculations, we need to 
determine the radial distribution of dose, from delta 
rays and the primary interactions, about the path 
of an energetic charged particle (Waligorski, Hamm, 
and Katz 1986). We presently use an analytic repre- 
sentation of the results of a Monte Carlo calculation 
made for liquid water for this purpose. More recently 
we have extended this model to include some solids 
used as radiation detectors (Katz et al. 1990). Ad- 
ditional information about both theoretical and ex- 
perimental determinations of the radial dose distri- 
bution may be found in Katz and Varma (1991). We 
use this information in connection with the response 
of the detector to gamma rays to find the radial dis- 
tribution of effect around a particle’s path. Because 
we interpret the response as the probability for acti- 
vating a target, we can make a map of the radial dis- 
tribution of activated targets. If, for example, we are 
interested in the opacity of a track in nuclear emul- 
sion, we can calculate the attenuation of a beam of 
light in a microscope photometer, as in the study of 
cosmic-ray tracks (Katz and Kobetich 1969). Alter- 
natively, we can integrate the probability radially to 
yield the cross section for the interaction of a single 
ion with the target (Katz 1978b; Waligorski. Loll, 
and Katz 1987). 

4.7 Hittedness 

In radiobiology we are not yet able to measure 
the effect produced by the interaction of a single ion 
with a target. Yet this ability to measure is of central 
importance in estimating the effects of GCR in space- 
flight. We can describe the interaction of a single ion 
with a target with the cross section. In track theory 
we wish to know whether a single particle, be it elec- 
tron, proton, a-particle, or whatever, is capable of 
inducing the tested end point with observable prob- 
ability. We characterize these interactions through 
the concept of hittedness , borrowed from biological 
target theory (Dertinger and Jung 1970). For a spe- 
cific irradiation the appropriate hittedness is either 
the number of interactions between charged parti- 
cles and target needed to induce the end point or 
the number of incident particles that must bombard 
the target, whichever is smaller. If either a single 
particle or a single interaction leads to the event, we 
will observe exponential response, as predicted by the 
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cumulative Poisson distribution. If two electrons are 
required for the event, we expect to observe a re- 
sponse described by the two- or- more- hit cumulative 
Poisson distribution. But the inactivation may take 
place through the transit of a single a-particle. In 
that case we expect to observe that the response to 
o-particles is one or more hit. We characterize the 
hittedness of a detector by its response to electrons 
or to gamma rays. Experimentally, if the response to 
the dose of gamma rays is exponential, we speak of 
a one-or-more-hit detector. 

4.8 The One-Hit Detector 

Most commonly, radiation detectors can be de- 
scribed as one-or-more-hit detectors. We imagine the 
detector to be a collection of targets- -sometimes ex- 
plicit, as in photographic emulsion, and sometimes 
implicit, as in a Fricke dosimeter or in alanine. Each 
of these targets is capable of responding to the tran- 
sit of a single electron of appropriate energy. The 
response is exponential; that is, it is linear at low 
dose and sublinear at high dose as the available tar- 
gets tend to have been inactivated. Saturation or 
overkill occurs at high dose. For one-hit detectors 
the response to heavy ions is also exponential with 
dose or fluence. 

To calculate the inactivation cross section for a 
one-hit detector, we first find P(D), the probability 
for target inactivation, as a function of the dose D of 
gamma rays. Next we fold this into the average ra- 
dial dose distribution about an ion’s path, to find the 
probability for target inactivation P(t) at radial dis- 
tance t. We integrate P(t) radially to find the action 
cross section a. When targets out to about three tar- 
get diameters are all inactivated, we simply use the 
point distribution of dose in our calculations. If we 
must take into account effects closer than three tar- 
get radii from the ion’s path it is necessary to average 
the dose in the extended targets to accommodate the 
dose gradient. Averaging is necessary because the ra- 
dial dose falls off essentially inversely as the square 
of the radial distance increases to a limit determined 
by the maximum delta- ray penetration. This limiting 
distance, determined essentially by the speed of the 
ion, places an upper limit on the action cross section 
that is observed experimentally with very heavy ions 
as thindown. ( Thindown is named after the appear- 
ance of the tracks of heavy ions in electron-sensitive 
emulsions, where the stopping end of a track looks 
like a sharpened pencil.) Typically for these detec- 
tors the cross section increases with an increase in 
(Z* / (5) 2 to a maximum (typically unrelated to tar- 
get size) and then declines in thindown. 


The RBE is equal to aE^/L, where Eq is the 1/e 
dose, or the dose for 37 percent survival. For one-hit 
detectors the RBE never exceeds one. The magni- 
tude of the cross section is approximately determined 
by the radial distance at which the dose equals Eq. 
We may speak of the cross section as approximating 
the size of the damaged region, but it is inappro- 
priate to speak of track size without specifying the 
end point. It is easy to estimate the cross section of a 
heavy ion with the grains of a nuclear emulsion from a 
microphotograph by estimating the radial distance at 
which about 63 percent of the grains are developed. 
At 63-percent development, we note that for insensi- 
tive emulsions where the track resembles a string of 
beads, the cross section is less than the grain size, 
whereas for a sensitive emulsion where the track re- 
sembles a hairy rope, the cross section may be orders 
of magnitude greater than the grain size. 

Our first venture into the one-hit detector was 
made for dry enzymes and viruses (Butts and Katz 

1967) . This work was followed by a model for the 
response of nuclear emulsions (Katz and Kobetich 
1969), of scintillation counters (Katz and Kobetich 

1968) , of thermoluminance detectors (TLD’s), of ala- 
nine (Waligorski et al. 1989), and for Escherichia 
coli B (Katz and Zachariah 1991). There are indi- 
cations that the plastic CR-39, used as an etchable 
track detector, is also a one-hit detector (Katz 1984). 

The global applicability of the model of the one- 
hit detector to a wide variety of detectors whose 
mechanisms arc vastly different from each other is 
at first thought to be rather astonishing. It arises 
simply from the fact that in each case the end point is 
stimulated by the passage of a single electron through 
the target volume. 

4.9 Supralinearity and the Linear 

Quadratic Model 

If a system has both one- hit and two- hit tar- 
gets having different radiosensitivities and popula- 
tions, we must expect that response will be linear at 
low dose, quadratic at intermediate doses, and sat- 
urating at high dose. We call the variable response 
supralinear and have proposed such a model to ex- 
plain supralinearity in TLD-100 (Katz 1978a). Note 
that the concept of a two- hit target requires only 
that two incident electrons are needed to stimulate 
the end point. The response may arise after pro- 
cessing as well as in the initial interactions, for the 
present model treats each detector as a black box. 
But if there are not two varieties of response, we 
cannot understand supralinearity in this model. Nor 
can we understand how a detector whose response to 
gamma rays is exponential can exhibit an RBE with 
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heavy ions greater than one unless there are tem- 
poral effects hinging on the time difference between 
gamma-ray and delta-ray exposures. Such a time 
difference between hits appears explicitly in a kinet- 
ics model that may provide an approach to temporal 
effects. 

In the same way, we do not understand the ra- 
tionale for using the linear quadratic formula to fit 
radiobiological data if there are not two types of tar- 
gets within a cell. We note that the formula is simply 
the first two terms of a series expansion and is usually 
applied to data of very limited dynamic range. Fur- 
ther we note that for radiations of different quality, 
there is no theory that can reliably predict the values 
of either the linear or the quadratic term. Neverthe- 
less, these terms are liberally interpreted with such 
phrases as “could be” or "might be,” though with 
equal validity one might insert “not.” 

Those who prefer to interpret data on the basis 
of hypothetical mechanisms whose details are rarely 
accessible via experiment may object to our paramet- 
ric formulations. Yet the track-structure model and 
its experimental parameters should not be dismissed 
lightly, for they may suggest mechanistic interpreta- 
tions that supersede those presently popular (Good- 
head 1989). We note that Newton’s laws were stimu- 
lated by Kepler’s phenomenology and that quantum 
theory was stimulated by Planck’s exercise in curve 
fitting. 

4.10 Cell Survival Model 

As detectors, biological cells require special con- 
sideration. For other detectors we assume, as in 
the case of nuclear emulsions, that the target has 
a characteristic size and no internal structure. The 
response of these detectors is then characterized by 
the parameters £q, the dose of gamma rays at which 
there is an average of one hit per target, «(), the sen- 
sitive site radius, and C y the hit ted ness. We also 
sometimes introduce a dimensionless track-structure 
parameter k proportional to E()Oy. Biological cells, 
however, have internal targets. We imagine the cells 
to resemble a bean bag in which the cell nucleus is the 
bag and the targets are the beans. We take it that 
the beans are one-hit in character but that m of the 
beans must be inactivated to generate the observed 
response. We also imagine that the beans are well 
distributed throughout the bag so that an energetic 
ion passing through the bag has the possibility for in- 
activating rn beans. Such a model makes it possible 
to understand why flatted cells respond differently to 
a-particles than rounded ones do. 

In this model the observed cross section is related 
to the size of the bean bag and the variation of 


response with LET is related to the properties of the 
beans. To set up a model of cellular response, we 
calculate the cross section for a hypothetical cluster 
of m overlapping beans and then assume that the 
cross section for the bean bag is proportional to 
that of the cluster (Katz, Sharma, and Homayoonfar 
1972). Because our model is based on the radial 
distribution of dose from delta rays, it automatically 
predicts thindown. For mammalian cells the fitted 
value of the k parameter suggests that the bean 
radius is about 1 /on, and thus the target for cell 
killing may be a chromosome. 

5. Mathematical Track-Structure Model 
of Cell Survival and Applications 

Detailed descriptions of the track model of cell 
survival have been given elsewhere (see, for exam- 
ple, Katz, Sharma, and Homayoonfar 1972). Here we 
present only the main concepts of the model and list 
the equations used in our calculations. Following our 
earlier studies of the appearance of particle tracks 
in nuclear emulsion (Katz and Kobetich 1969), the 
model distinguishes between the grain-count regime, 
where inactivations occur randomly along the par- 
ticle’s path, and the track-width regime, where the 
inactivations are distributed like a hairy rope. The 
transition from the grain-count to track- width regime 
takes place at Z* 2 /k,/3 2 ss 4; lower values are in the 
grain-count regime and higher values arc in the track- 
width regime. The quality k is a parameter of the 
track-structure model that combines both the tar- 
get size and the characteristic dose of gamma rays at 
which there is an average of one hit per target. As 
in nuclear emulsions we speak of a thindown regime 
where the cross section is limited by the kinematic 
constraint on delta-ray energies, but has nothing to 
do with the Bragg peak in stopping power or with 
the changing effective charge of a slowing-down ion. 

To accommodate for the capacity of cells to ac- 
cumulate sublethal damage, two modes of inactiva- 
tion are identified, namely ion kill (intratrack) and 
gamma kill (intertrack). In these two inactivation 
modes, the statistical character of the inactivation 
changes rather than the fundamental physical inter- 
action. Effects are in terms of dose rather than the 
number of electrons passing through the nucleus. We 
do not find justification for considering the stopping 
end of an electron track as a source of ion kill nor do 
we consider the radial separation of a heavy-ion t rack 
into core producing ion kill and penumbra producing 
gamma kill. 

The cell survival model uses Z* 2 //? 2 as a plotting 
parameter superior to LET, now in wide use. At the 
stopping end of a track, at highest LET, even this 
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parameter fails because in the thindown regime the 
cross section depends on /?, the relative speed of the 
ion. In the thindown regime the cross section is some- 
times plotted against energy per unit mass, a related 
parameter. The model bases the meaning of low LET 
on the comparison of Z* 2 /,# 2 with k. Similarly, low 
dose means low compared with E[). The model ex- 
plains why plots of extrapolated cross section from 
the tail of a survival curve tend to be single-valued 
functions of LET at low LET (response is dominated 
by gamma kill) and why they are multiple valued 
with Z at high LET (response is dominated by ion 
kill and thindown). It explains why plots of RBE 
versus LET for biological cells pass through a maxi- 
mum when about half the intersected cells are killed 
in ion kill. It predicts that the RBE for lighter ions 
will be greater than the RBE for heavy ions at the 
same LET and the same survival level. This inequal- 
ity is because of the structure of particle tracks. At 
the same LET, the heavier ions move faster. Its delta 
rays are fewer but more energetic and gamma kill is 
more likely, which reduces the RBE. 


The gamma-kill dose is 


D 7 - (1 - P)D 


where the ion-kill probability is 



(4) 


In the track-width regime, where P > 0.98, we 
take 

II, = 1 (0) 

and find a from the track width, which increases 
linearly with Z* f A as the inactivation cross section 
increases with ZD 2 j A 2 up to the limit set by the 
maximum radial range of delta rays. This is the 
thindown region. 

To find the cross section in the track-width 
regime, including the thindown region, a. separate 
calculation must, be made. First the target cross sec- 
tion S for targets of radius uq is found from 


5.1 Gamma Kill 

Cells not inactivated in the ion-kill mode can be 
sublethally damaged by delta rays from the passing 
particle and then inactivated in the gamma-kill mode 
by cumulative addition of sublethal damage from 
delta rays from other passing ions. Survival in the 
gamma-kill mode is taken to follow the m- target 
statistics of inactivation by secondary electrons from 
X-ray or gamma- ray photons. 

5.2 Mathematical Formalism 

In the grain-count, regime the surviving fraction of 
a cellular population whose radiosensitivity parame- 
ters are rn, (7 q, and k, after track-segment irradi- 
ation with an ion dose D of a fiuence of F particles of 
charge number Z, effective charge value Z*, relative 
speed /3j and stopping power L (LET :X; ), is found 
from the expression 


N 

Ao 


IL x IT 


where ion-kill fraction is 

II,: = exp(— aF) 
and where gamma-kill fraction is 

-Eh 


n 7 = i 


1 — exp 


Eo 


(i) 


(2) 


(3) 


S — 2n 


P l <lt ( l - 
Jo 


where D(t , ao) is the average' radial dose in the target 
of radius «(j. The target radius ao is found from k 
and E{) according to 


k — (E()d 2 ) x (fix lO^erg-em) 


(8) 


with multitarget response to gamma rays character- 
ized by E{) and m found for the cell. The S value 
must be multiplied by the ratio of the plateau value 
of the action cross section to the plateau value of 
the target cross section S{\ to yield the cellular action 
cross section in the track- width regime (Katz et ah 
1971). In this region we assume there is no gainma- 
kill dose, though in the outer reaches of the track 
width some small fraction of the energy lost, by the 
ion is deposited in the gamma-kill mode. 

To calculate RBE at a given kill (transformation) 
or survival level we use the definition 


RBE = T 


where 


D r = — Eo < In 


N \ C”' 
No 


(9) 


( 10 ) 


is the equivalent X-ray dose, and D is the correspond- 
ing ion dose. 
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All calculations pertain to water, so the ion dose 
is always 

D = FL (11) 

Where cross sections and RBE’s are calculated from 
the final slope of the survival curves, we refer to the 
cross section and RBE as extrapolated, and in the 
grain-count regime we write 


a ext — a oP + ^ (12) 

and 

RBE ext = ( * q Eq/L)P + (1 -P) (13) 

To calculate the effective charge value of an ion of 
charge number Z moving with a relative velocity /?, 
we use the expression (Barkas 1963) 

Z* = Z[ 1 - exp(-125 (3Z~ 2/3 )\ (14) 

We calculate the stopping power and range in water 
of an ion of atomic number Z with the expression 

L(Z,E) = L(p,E) (E\ (15) 


beams (Roth and Katz 1980), and most recently for 
cosmic rays (Cucinotta et al. 1990). 

5.3 Cell Killing, Chromosome 

Aberrations, and Transformations 

Our treatment of transformations (Waligorski, 
Sinclair, and Katz 1987) is based on data obtained 
with the BEVALAC accelerator by Yang et al. 
(1985). We use the same form of equations for cell 
killing and take cell killing and transformation to be 
independent processes that take place along the same 
particle track. Our parameters for Chinese ham- 
ster cells are based on the data of Skarsgard et al. 
(1967) whereas the parameters for Tradescantia are 
based on the data of Underbrink et al. (1978). We 
have extracted parameters from data obtained with 
X rays and neutron irradiations of two different en- 
ergy spectra, ignoring possible gamma-ray contami- 
nation. The parameters for T-l cells of human ori- 
gin are taken from Todd (1967). The separate sets 
of parameters for cell survival, chromosome aberra- 
tion, and transformation are shown in table 1 (Katz 
and Huang 1991). Two sets of parameters are shown 
when the data do not permit a clear distinction be- 
tween them. 


where Z* and Z* are the effective charges of the ion 
and proton, respectively, and L(p , E) is the stopping 
power, in water, of a proton at the same energy per 
nucleon E. 

At low fluence, where ions are sufficiently far 
apart that intertrack effects are unlikely, we can 
neglect the contribution from gamma kill. Under this 
circumstance RBE is 

RBE = E 0 (^) 1/m (16) 

This value is applicable to low doses of neutrons as 
well as to the effects of GCR’s (Katz and Cucinotta 
1991) and agrees with the measurements of Bettega 
et al. (1990) down to doses of 0.01 Gy. 

In order to determine the RBE for mixed radi- 
ation fields, our model requires knowledge of the 
particle-energy spectrum of the radiation field. We 
then find the totality of effects due to ion kill, add 
to that sum the gamma-kill doses including the dose 
from gamma rays, find from this value the ion-kill 
and the gamma-kill survival probability, and take 
their product to be the surviving fraction of irradi- 
ated cells. We have done these calculations for neu- 
trons admixed with gamma rays (Katz, Sharma, and 
Homayoonfar 1972), for range-modulated heavy-ion 


5.4 Target Fragmentation Effects 

High-energy protons passing through tissue will 
occasionally cause nuclear reactions that produce 
low-energy, high-LET ions from the tissue itself. The 
target fragments, in turn, will be a source of delta 
rays that should contribute to biological damage lo- 
cally in the tissue matrix. The differential fluence Fj 
(Wilson 1977) describes the local source of target 
fragments: 


J L(Zj,E) 



dEj(E') 

dE 1 


Fp(Ep) dE' 


(17) 


where j is the fragment label, L is the stopping 
power, Ej is the macroscopic nuclear production 
cross section, and F p is the fluence of protons with 
energy E p . An effective action cross section for the 
proton dressed by the local target fragments (nuclear 
stars) is now written as 


2 roc 

a ~ a p {E p ) + F /p \ XI / Fj(Ej) crj(Ej) dEj 

r p\^p) ■ j 0 

(18) 

where a p and Oj are given by the Katz formalism 
(eq. (5)). The gamma-kill dose for the proton and 
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target fragments is written as 

D 1 = D^pit'p) + ^ f Fj(Ej) L(Zj, Ej) [l - PjiEj)] dE } 
j J " 

(19) 

The production energy spectra for the target frag- 
ments are expressed as (Wilson et al. 1989) 


dE 


Ej(.Ep) \fE E / 2 £ 0j 

Mf 2 


( 20 ) 


where the average energy of the fragment is given 
by 3E 0j . 

The fragmentation parameters that arc described 
in this paper are discussed in Wilson, Townsend, and 
Khan (1989). The light-ion production cross sections 
are from Bertinrs Monte Carlo results (Anon. 1968) 
and the Silberberg-Tsao empirical model is used for 
the heavier fragments (Silberberg, Tsao, and Shapiro 
1976). The average energy of the tissue fragments 
is related to the momentum width that is measured 
experimentally (Greiner et ah 1975), which Wilson, 
Townsend, and Khan (1989) fit empirically. The 
largest uncertainties exist for light-ion production 
(A = 7 and 9) and for energies below 100 MeV. 
Elastic recoils and meson production above several 
hundred MeV will also contribute to the action cross 
section and should be added. Stopping power in 
tissue is from the work of Wilson (1983), which is 
based on the Ziegler (1980) analysis. 

The solid line in figure 1 displays the frag- 
ment LET component derived from equation (17) for 
1-GeV protons in water. The dotted and dashed lines 
show the contributions from the secondary proton 
fragments and a- particles, respectively. 

The proton action cross section for cell survival of 
V79 Chinese hamster cells versus the proton energy is 
shown in figure 2. The cellular response parameters 
are given in table 1. These parameters are compared 
with the data of Hall et al. (1978) for survival of 
Chinese hamster cells where the characteristic X-ray 
dose Eq is taken as 2.9 Gy, as found from their X-ray 
data and with the remaining parameters the same as 
given by table 1. The dotted line in figure 2 shows the 
contributions from primary ionizations; the dashed 
line shows contributions from secondary ions; and 
the solid line shows total ion contributions. The 
oxygen and nitrogen fragments contribute partially 
to the cross section in the track-width regime. The 
decreasing proton LET with increasing energy leads 
to complete domination by target fragments above 
about 50 MeV. The shape of the action cross section 


in figure 2 directly reflects the nuclear absorption 
cross section in tissue. We expect a further increase 
above several hundred MeV when meson production 
is included in the cross section. In table 2, the 
individual contributions to the action cross section 
for V79 Chinese hamster cells are shown for several 
proton energies. Secondary protons and o-particles 
arc dominant, with broad spectra of tissue fragments 
making nonnegligible contributions. Primary proton 
contributions make up an insignificant fraction of 
the action cross section above 100 MeV, and the 
relatively slow change with energy of the nuclear 
production cross sections leads to a plateau in the 
action cross section at high proton energies. 

Figure 3 shows the action cross section versus 
proton LET with the calculations of figure 2 corre- 
sponding to high-LET protons depicted to 0.1 MeV. 
The behavior of the cross section below 0.5 keV//iin 
shows the dominance of the tissue secondaries (nu- 
clear stars). At about 0.2 keV//xm, the proton LET 
reaches a minimum and then increases, which is the 
origin of the ‘hook” in figure 3 at the lowest LET 
values. It would be interesting to test our results 
for the proton cross section by experiment. Results 
herein assume an equilibrium in the local secondary 
fluence spectra and are sensitive to interface effects 
(Cucinotta, Hajnal, and Wilson 1990) and the com- 
position of the host media of the cell culture. 

Figure 4 shows the proton gamma-kill dose di- 
vided by Eq versus proton energy. The primary ion- 
ization is the dotted line and the solid line includes 
the effects of fragments. Secondary ion production 
has a negligible effect on intertrack effects, except at 
the highest energies, where a small contribution is 
seen. 

Cellular parameters obtained for survival and 
neoplastic transformations of C3H10TV2 cells ob- 
tained from the experiments of Yang et al. (1985) 
are given in table 1. The large uncertainties in the 
transformation data of Yang et al. in table 1 should 
lead to similar uncertainties in the transformation 
parameters. Parameter sets were found from data for 
instantaneous and delayed plating of the cells after 
irradiation. In this paper only the delayed plating 
case is considered. Using the parameter sets, gen- 
eral agreement with the measured RBE values was 
found (Waligorski, Sinclair, and Katz 1987). The 
single-particle-inactivation cross section with the tar- 
get fragmentation of equation (18) neglected is shown 
in figures 5 and 6 for cell death and cell transforma- 
tion, respectively, as a function of the energy of the 
passing ion. The target fragmentation contribution 
(the second term in eq. (18)) for protons has been 
evaluated as shown in figures 7 and 8. For protons 
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the effects of the target fragments, shown as a dashed 
line, dominate over the proton direct ionization (dot- 
ted line) at high energy. For high-LET particles (low 
energy), direct ionization dominates and target frag- 
mentation effects become negligible. A simple scaling 

by 2 relates the proton target fragment term to 
ions of mass Aj. The resulting effective action cross 
sections for cell death and cell transformation are 
plotted in figures 9 and 10, respectively. We note that 
the low-energy J ^Fe component of the GCR spectra 
extends into the track- width regime where o > 
and is not represented in the present calculations. 
The resultant error introduced into the present cal- 
culation is small. 

5.5 Survival Curves and Proton RBE 

Cell survival curves for 10-, 100-, and 1000-MeV 
protons are shown as a function of absorbed dose 
in figures ll(a)-(c). The dotted lines shown in the 
figures indicate the primary ionization contributions 
and the solid lines indicate the added effects of the 
target fragment contributions. Results show the im- 
portance of secondary ion production for increasing 
energy. For example, at 1000 MeV the increase in 
cell death due to the fragments docs not lead to sub- 
stantial changes in RBE at high doses, as shown in 
figure 11(c). It is in the initial portion of the sur- 
vival curves whore the ion-kill mode causes large dif- 
ferences in RBE when compared with gamma rays. 

Proton RBE for survival of Chinese hamster cells 
versus absorbed dose is shown in figure 12, and all 
curves include the effects of target fragmentation. 
We note that in figure 12, proton fluence is F p = 
6.24D/LET, with F p in protons/^m 2 . The rise in 
RBE at low dose or fluence, where single proton 
tracks dominate, is directly attributed to ion kill 
from both primary protons at 10 MeV and nuclear 
fragments at the higher energies. Not shown are RBE 
calculations neglecting the target fragments, which 
are nearly identical to the results for the 10-MeV 
proton in figure 12 and arc almost identical to unity 
for the 100- and 1000-MeV protons. The low-dose 
behavior of the RBE can be seen from equation (16), 
where a dependence on Z) -2 ^ is found for m = 3 
(from the data of Skarsgard et al. (1967) in table 1). 
This effect is supported experimentally as discussed 
below (Cucinotta et al. 1991c). 

5.6 RBE of 160-MeV Plateau Region 

Protons 

Cell survival experiments have been performed at 
the Harvard Cyclotron for the purpose of determining 
the RBE of the protons. The V79 Chinese hamster 


cells cultured in vitro were irradiated in the plateau 
region of the Bragg curve and in a spread-out Bragg 
peak by Hall et al. (1978). Using the track model and 
the high-energy nucleon transport code BRYNTRN 
of Wilson et al. (1988 and 1989), we compare the 
survival measurements and RBE determinations for 
attached cells in the plateau region of the 160-MeV- 
proton Bragg curve. 

The nucleon transport code BRYNTRN solves 
the coupled proton-neutron transport problem for 
high energies in the straight-ahead approximation 
with multiple-scattering and straggling effects ig- 
nored. Target fragments with A > 1 are transported 
by taking the production collision density, which is 
given as Fj in equation (17). The respective Boltz- 
mann equations for proton and neutron transport are 




$ p (x, E ) 

roo 

f pj {E,E')* j (x,E')dE , { 21) 


and 


d 

dx 


+ Z n (E)\<i>n(x,E) 


y-.i: f n j ( E, /•;') *j{x,Ef) dE' (22) 


where )j is the particle flux of type j particles at po- 
sition x with energy E\ L(Zj , E) is the proton stop- 
ping power; T^ p {E) and E n (E) are proton and neu- 
tron total cross sections, respectively; and fij{E, E ! ) 
represents the differential cross sections for elastic 
and inelastic processes. As described by Wilson 
and Lamkin (1975), the Boltzmann equations (21) 
and (22) are solved with a characteristic transforma- 
tion to reduce the problem to a set of coupled inte- 
gral equations with boundary conditions at x = 0, 
which are then solved numerically. More details on 
the method of solution and the nuclear scattering 
data base are given by Wilson et al. (1988 and 1989). 
The Bragg curve obtained from BRYNTRN for 160- 
MeV protons in water is shown in figure 13 with mea- 
surements of Verhey et al. (1979). In figure 3 3 the 
squares represent the primary dose and the circles the 
total dose including secondary production. Calcula- 
tions are normalized to the peak of the experimental 
Bragg curve. Straggling and multiple-scattering ef- 
fects, which are not included here, both contribute 
significantly at the peak of the Bragg curve. We 
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consider the plateau region to be where the high- 
energy assumptions are approximately true. 

At energies of 160 MeV, nuclear recoils from elas- 
tic scattering provide a sizable correction to the sec- 
ondary ion production represented by target frag- 
mentation. Elastic nuclear scattering is represented 
by the Born term of the optical model renormal- 
ized to the total scattering cross section in the 
BRYNTRN code. This representation of elastic nu- 
clear scattering is fairly accurate for integral quanti- 
ties above 100 MeV, but it breaks down at lower en- 
ergies because of multiple scattering, nuclear medium 
corrections, and especially Coulomb effects. The cor- 
rection to the proton action cross section from clastic 
scattering is shown in table 3 for several energies and 
is included in the following comparisons. 

Results for the surviving fraction of suspended 
V79 Chinese hamster cells irradiated by 60 Co 
gamma rays and plateau region of a 160- MeV proton 
beam (Hall et al. 1978) are shown in figure 14. The 
dashed line is the fit to the experimental gamma-ray 
survival curve, the dotted line (barely distinguish- 
able from the gamma- ray response) is the contribu- 
tion from primary protons only, and the solid line 
depicts calculations that include the effects of nu- 
clear reactions. The characteristic gamma-ray dose 
E{ ) is taken as 2.9 Gy to reproduce the experimental 
gamma-ray curve with the other response parame- 
ters given above. The dashed and dotted lines are 
nearly identical, an indication that high-energy pro- 
tons minus the effects of nuclear force indeed act as 
gamma rays. Agreement with the data is fair and 
indicates that the modeling of nuclear fragmentation 
made in this paper is somewhat lacking. Figure 15 
presents our results for the proton RBE and the val- 
ues obtained with the analysis methods of Kellerer 
and Brenot (1973) as discussed by Hall et al. (1978). 
The “bare” proton RBE has a value of 1 (not shown), 
except at the lowest doses, where a small contribution 
from ion kill gives a slight increase. The calculations 
of the proton RBE presented in figure 15 show good 
agreement with the experimentally obtained values, 
with the increasing RBE at low doses shifted to lower 
values than experiment. A second analysis method- 
ology that assumes only that the dose-effect, curve 
is convex but is otherwise shape independent was 
used to derive RBE values (Hall et al. 1978) and they 
are shown in figure 16 with present predictions. The 
rise in RBE values at low dose as predicted in equa- 
tion (16) is clearly seen in the calculations and the 
experimental values. Also shown in figure 16 are our 
calculations neglecting nuclear reactions, the values 
being almost exactly 1 for all doses. 


5.7 Cell Survival in HZE Beams 

The HZE transport problem has been solved and 
related to the Bragg curve (Wilson 1977 and 1983) for 
monoenergetic unidirectional ion beams. The Bragg 
curves we calculate also provide the values for fiuenee 
estimates for the exposure conditions of biological 
samples yet to be analyzed. Errors in the Bragg 
curve translate directly into errors in exposure levels 
for comparisons with experimental response data. 

Calculations of the relative ionization ratio for 
^Ar at 514 MeV/nucleon are compared with the ex- 
periments of Blakely et al. (1979) in figure 17. The 
effects of nuclear secondaries arc 1 seen to be most 
important near and beyond the Bragg peak. Sur- 
vival for aerobic and hypoxic T-l cells of human ori- 
gin has been calculated with the Katz parameters 
in table 1 for several locations along the beam line 
within a water column for three different ion beams 
of C, Ne, and Ar. Calculations include both projec- 
tile as well as target fragments. Results for the 
beam experiments (Blakely et al. 1979) are shown 
in figure 18. Effects of overlapping delta rays are 
clearly apparent except near the Bragg peak, where 
the sigmoid appearance has all but disappeared. The 
sigmoid shape returns downstream from the Bragg 
peak, where the overlapping delta-rays from adja- 
cent ions again contribute to the exposure. Note 
that the oxygen effect has all but vanished near the 
Bragg peak whereas hypoxic cells show considerable 
radiation resistance both upstream and downstream 
from the Bragg peak. Results of our calculations for 
2 ^Ne beams are shown in figure 19. Results for the 
2() Ne beam arc qualitatively similar to those for 
beams. The region over which t he sigmoid appear- 
ance is suppressed is greatly expanded in the Bragg 
peak region. Oxygen enhancement is greatly dimin- 
ished 1 cm before and after the Bragg peak, as can 
be seen in figure 20. This fact is of potential impor- 
tance to radiation therapy. The sigmoid behavior is 
virtually nonexistent for the 4t) Ar beam exposures, as 
shown in figure 20. Obviously, at some great distance 
downstream, the sigmoid shape will appear because 
only light fragments will survive. No experiments 
were conducted in this region. 

5.8 Cell Damage for the GCR Spectrum 

To apply the cellular track model to the mixed- 
radiation fields found in space, we not'd to make 
the appropriate replacement of the cross section and 
particle fiuenee number aF with the particle field 
quantities and their corresponding cross sections. 
The ion-kill term, which will now contain a projectile 
source term (including projectile fragments) and a 
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target fragment term, is written as 

aF = Y, [ dE l ^ x <Ej) 

3 J 

+ EL / dE n dEj $ a (x, E„ : Ej) a n (E n ) (23) 

a j 

where the second term is the contribution of nuclear 
fragments produced locally in the biological medium 
(Cucinotta et al. 1991a). This contribution may 
also be written in terms of an effective action cross 
section cr* for the passing ion, whose track is dressed 
by the local target fragments (nuclear stars), as 

oF = £ J dE i ( 24 ) 

j 

The gamma-kill dose fraction becomes 

D 7 = Z / dE i *A x > E j) I 1 - p j( E j)) s A e j) 

j 

+ H [ dE j dE a *a(®, Ea : Ej) 
j 

x [1 -Pa(Ea)} Sa(E a ) (25) 

Equations (23) and (25) are used in equations (2) 
and (3), respectively. The summations over all par- 
ticle types in equations (23) and (25) represent the 
addition of probabilities from all ions in the radiation 
field that contribute to the end point under study. 

The cellular track model was applied to predict 
the fraction of C3H10TV2 cells killed or transformed 
for 1 yr in deep space at solar minimum behind 
typical spacecraft shielding. The GCR environment 
was taken from the Naval Research Laboratory code 
(Adams, Silberberg, and Tsao 1981). Aluminum 
shielding was considered with a local region of tis- 
sue for the cell cultures. Tables 4 and 5 contain 
individual particle fluences and absorbed doses, re- 
spectively, for the protons, a-particles, Z = 3 to 
9 ions (labeled L-Z), and Z = 10 to 28 ions (labeled 
H-Z) as determined by the Langley GCR code 
(Wilson, Townsend, and Badavi 1987). Results for 
the fraction of C3H10TV2 cells killed and transformed 
for 1 yr at solar minimum behind aluminum shield- 
ing are listed in tables 6 and 7, respectively. The 
gamma-kill mode was of negligible importance in the 
calculations, and this unimportance indicates that 
biological damage in deep space from GCR particles 
at the cellular level will indeed result from the ac- 
tion of single particles. The importance of the target 


terms in biological effects for low-LET protons and a- 
particles is quite apparent. The results also indicate 
that the HZE component of the GCR spectrum is the 
most damaging for small shielding depths. At large 
depths the HZE components break up and cause pro- 
ton buildup with increasing shield depth. At large 
depths the protons dominate the biological effects. 
In comparing individual charge components, we see 
that the H-Z particles have a reduced effectiveness 
for the transformation end point. 

Also listed in tables 6 and 7 are the values of RBE 
for the two end points. In table 8 we show the present 
RBE values beside the average QF values taken from 
Townsend et al. (1990), who used the same transport 
code. The fact that RBE and QF values are nearly 
equal at small depths is somewhat coincidental. We 
note that the QF is independent of the fluence level; 
this independence is not true for the Katz model. 
The Katz model indicates a substantial increase in 
risk, at higher shielding levels, compared with QF’s 
given in report no. 26 of the International Committee 
on Radiological Protection (Anon. 1977). 

The RBE values show a simple scaling with expo- 
sure time for the GCR particles, as can be seen from 
equations (9), (10), and (2) when ion kill dominates. 
Here we find for 



N 0 

(26) 

with 


oF <$: l 

(27) 

that 


RBE = ^1/™ p'l-i+ii/™')] 

L 

(28) 


Then, scaling RBE as a function of duration in deep 
space to the 1-yr value r\ for a duration period of 
r (with F — nr) gives (Cucinotta et al. 1991a and 
1991b) 


/ \ [- 1 + 0 /™)] 

RBE(r)=^-J RBE(ri) (29) 


As a result, a one-hit (m =1) system RBE becomes 
fluence independent, as expressed by 

RBE(r) = RBE(n) (30) 

a two-hit (m = 2) system is expressed by 


RBE(r) 


RBE(n) 

(t/ti ) 1 / 2 


(31) 
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and a three-hit (m = 3) system is expressed by 


RBE(t) 


RBE(ti) 

(r/n) 2 / 3 


(32) 


Results of this scaling approximation agree quite well 
with calculations from the Katz model, as shown 
in table 9, where values obtained with the approx- 
imations of equation (28) are shown in parentheses 
as scaled from the 1-yr RBE values taken from ta- 
ble 8, and results of the calculations are shown with- 
out parentheses. The extremely large RBE values 
that would be obtained for small values of r are due 
to the choice of energetic photons as the reference 
radiation. 


6. Concluding Remarks 

Over the past 25 years Katz and coworkers have 
developed a model of particle tracks that began with 
nuclear emulsions and subsequently has been ex- 
tended to other detectors and to the biological ef- 
fects of high linear energy transfer (LET) radiations. 
This model requires as input information knowledge 
of the particle-energy spectrum of the radiation envi- 
ronment as well as the dose of gamma rays. Calcula- 
tions of the effects of beams of protons, of heavy ions, 
and of energetic neutrons have been hindered because 
of the lack of a model of such beams that included 
both projectile and target fragmentation. A beam 
model created at NASA Langley Research Center 
remedied this neglect. Through it we have been able 
to validate both the track theory of biological effects 
and the beam model by comparison of our calculated 
radiobiological end points with ground-based mea- 
surements for proton and heavy-ion beams. Based 
on this validation we have initiated calculations of 
biological effects in space vehicles in selected orbits, 
incorporating knowledge of the distribution of solar 
and galactic cosmic rays to be encountered there. We 
know of no other way to estimate the biological dam- 
age in space flight at the very low fluences of heavy 
ions to be encountered there. 


NASA Langley Research Center 

Hampton, VA 23681-0001 

August 20, 1992 
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Table 1. Parameters for Cell Survival, Aberrations, and Transformations 


Biological end point 

rn 

K 

Eq, 

Gy 

9 

cm 

CH 2 B 2 Chinese hamster cells 
(Skarsgard et al. 1967): 

Survival 

3 

1100 

1.82 

4.3 x 10" 7 

Abnormal metaphases .... 

3 

900 

1.82 

3.0 x 10" 7 

Chromatid exchanges 

2 

1400 

25 

6.5 x 10~ 9 

C3H10TV2 mouse cells 
(Yang et al. 1985): 

Survival 

3 

750 

1.7 

5.0 x 10" 7 

Transformations 

2 

750 

180 

1.2 x 10~ 1Q 

Transformations 

3 

475 

50 

7.0 x 10" 11 

Tradescantia 

(Underbrink et al. 1978): 

Survival 

2 

1000 

2.1 

3.5 x 10" 7 

Survival 

1.5 

1900 

2.6 

4.0 x 10" 7 

Human T-l cells 
(Todd 1967): 

Survival (aerobic) 

2.5 

1000 

1.7 

6.7 x 10" 7 

Survival (hypoxic) 

2.5 

1300 (1450)“ 

4.6 (5.2)“ 

6.7 x 10" 7 


a Two sets of parameters are shown because the data do not permit a clear distinction between 
them. 




Table 2. Target Fragment Contributions to Proton Action Cross Section for 

V79 Chinese Hamster Cells 





Table 3. Secondary Ion Production Contribution to Proton Action 
Cross Section for V79 Chinese Hamster Cells 


E p , MeV 

Contribution to proton action cross section. <7, 
cm 2 , for- 

Elastic recoils 

Fragmentation 

100 

1.13 x 10” 11 

1.081 x 10" u 

150 

.82 

1.02 

1000 

.58 

2.09 


Table 4. Fluence for 1 Year at Solar Minimum Behind Aluminum Shielding 


x , g/cm 2 

Fluence, particles/cm 2 , from 



Protons 

Q-particles 

"L-Z 

6 H-Z 

0 

1.29 x 10 8 

1.24 x 10 7 

1.09 x 10 7 

3.0 x 10 7 

1 

1.31 

1.21 

1.05 

2.8 

2 

1.33 

1.18 

1.01 

2.7 

3 

1.34 

1.15 

.98 

2.5 

5 

1.36 

1.10 

.91 

2.2 

10 

1.40 

.97 

.77 

1.7 

20 

1.43 

.77 

.57 

1.1 


a Z — 3 to 9 ions. 
b Z = 10 to 28 ions. 


Table 5. Absorbed Dose for 1 Year at Solar Minimum Behind Aluminum Shielding 


x , g/cm 2 


Dose, cGy/yr, from 

— 


Protons 

a-particles 

“ L-Z 

b H-Z 

Total 

0 

6.2 

3.0 

2.8 

5.0 

17.1 

i 

6.3 

2.7 

2.5 

3.6 

15.1 

2 

6.8 

2.6 

2.4 

3.3 

15.1 

3 

7.1 

2.6 

2.3 

3.1 

15.0 

5 

7.6 

2.4 

2.1 

2.7 

14.8 

10 

8.5 

2.1 

1.7 

2.0 

14.3 

20 

9.5 

1.7 

1.1 

1.1 

13.4 


a Z = 3 to 9 ions. 
b Z = 10 to 28 ions. 
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Table 6. Fraction of C3H10TV 2 Cells Killed in Deep Space for 1 Year at Solar Minimum 

Behind Aluminum Shielding 


x,g/cm 2 

Fraction of cells killed by — 

Total 

RBE 

Protons 

a-par tides 

° L-Z 

b H-Z 


Including target fragments 


0 

1.35 x 10~ 2 

0.46 x 10 -2 

0.57 x 10~ 2 

2.08 x 10 2 

4.46 x 10 -2 

7.1 

1 

.76 

.15 

.43 

1.84 

3.18 

7.0 

2 

.80 

.14 

.41 

1.69 

3.04 

6.9 

3 

.83 

.14 

.38 

1.55 

2.90 

6.8 

5 

.88 

.14 

.34 

1.32 

2.68 

6.7 

10 

.95 

.12 

.25 

.91 

2.22 

6.5 

20 

1.02 

.09 

.15 

.49 

1.74 

6.2 

Without target fragments 


0 

0.84 x 10~ 2 

0.37 x 10 -2 

0.55 x 10“ 2 

2.08 x 10“ 2 

3.79 x 10 -2 

6.7 

1 

.24 

.06 

.41 

1.83 

2.54 

6.5 

2 

.28 

.06 

.39 

1.68 

2.41 

6.3 

3 

.31 

.06 

.37 

1.55 

2.27 

6.2 

5 

.35 

.06 

.33 

1.31 

2.04 

6.1 

10 

.42 

.05 

.24 

.91 

1.61 

5.7 

20 

.49 

.04 

.14 

.48 

1.15 

5.3 


a Z = 3 to 9 ions. 
b Z = 10 to 28 ions. 


Table 7. Fraction of C3H10TV2 Cells Transformed in Deep Space for 1 Year at Solar Minimum 

Behind Aluminum Shielding 


g/cm 2 

Fraction of cells transformed — 

Total 

RBE 

Protons 

a-particles 

a L-Z 

b H-Z 

Including target, fragments 

0 

5.2 x 10~ 6 

2.0 x 10 -6 

3.1 x 10“ 6 

7.5 x 10“ 6 

1.78 x 10 -5 

6.4 

1 

3.5 

1.0 

2.7 

6.7 

1.39 

6.4 

2 

3.7 

1.0 

2.6 

6.2 

1.35 

6.3 

3 

3.9 

.9 

2.4 

5.7 

1.29 

6.3 

5 

4.2 

.9 

2.2 

4.9 

1.22 

6.2 

10 

4.7 

.8 

1.7 

3.5 

1.06 

6.0 

20 

5.2 

.6 

1.1 

2.0 

.88 

5.7 

Without target fragments 


0 

3.2 x 10 -6 

1.6 x 10^ 6 

3.1 x 10~ 6 

7.5 x 10“ 6 

1.53 x 10- 5 

6.0 

1 

1.4 

.6 

2.7 

6.7 

1.13 

5.8 

2 

1.6 

.6 

2.5 

6.2 

1.09 

5.7 

3 

1.8 

.6 

2.4 

5.7 

1.05 

5.6 

5 

2.1 

.5 

2.1 

4.9 

.97 

5.4 

10 

2.5 

.5 

1.6 

3.5 

.82 

5.2 

20 

3.0 

.4 

1.0 

2.0 

.64 

4.9 


a Z = 3 to 9 ions. 
h Z = 10 to 28 ions. 
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Table 8. Comparison of Average Quality Factors 0 Versus RBE for 
Cell Death and Transformation 

[One year in deep space at solar minimum] 


x, g/cm 2 

QF 

(b) 

RBE for cell 
death 

RBE for cell 
transformation 

0 

7.1 

7.1 

6.4 

1 

5.6 

7.0 

6.4 

2 

5.3 

6.9 

6.3 

3 

5.1 

6.8 

6.3 

5 

4.7 

6.7 

6.2 

10 

3.9 

6.5 

6.0 

20 

3.2 

6.2 

5.7 


0 Townsend et al. (1990). 
6 ICRP 26 (Anon. 1977). 


Table 9. RBE for Cell Death and Transformation of C3H10TV 2 Cell for GCR 
Spectrum at Solar Minimum Behind Aluminum Shielding 

[Values in parentheses scaled from 1-year value using cq. (29)] 


x , g/cm 2 

RBE values for time periods of 

1 mo 

1 yr 

2 yr 

Cell death 

0 

33.2 (37.0) 

7.1 

4.8 (4.6) 

1 

33.2 (36.1) 

7.0 

4.7 (4.5) 

3 

32.4 (35.1) 

6.8 

4.5 (4.3) 

Cell transformation 

0 

22.3 (22.2) 

6.4 

4.6 (4.5) 

1 

22.0 (22.2) 

6.4 

4.5 (4.5) 

3 

21.6 (21.8) 

6.3 

4.4 (4.4) 
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Figure 1. Integral LET spectra derived from equation (17) 
for nuclear fragments produced by 1-GeV protons in 
water. 


Total contributions 

Secondary contributions 



Figure 3. Calculated values of the proton action cross section 
for survival of V79 Chinese hamster cells versus proton 
LET. 



Figure 2. Calculated values of the proton action cross section 
for survival of V79 Chinese hamster cells versus proton 
energy. 


Total contributions 
(includes fragmentation) 
Primary ionization contributions 



Figure 4. Calculated values of the proton gamma-kill dose 
for survival of V79 Chinese hamster cells versus proton 
energy. 
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Figure 5. Cell-death action cross section for various ions 
in C3H10TV2 cells according to Katz model for direct 
ionization effects only. 



Figure 7. Cell-death action cross section including effects of 
nuclear reactions for protons in C3H10TV2 cells according 
to Katz model. 



Ep, MeV/nucleon 

Figure 6. Cell-transformation action cross section for various 
ions in C 3 H 10 TV 2 cells according to Katz model. 



Figure 8. Cell- transformation action cross section inc luding 
effects of nuclear reactions for protons in C3H10TV2 cells 
according to Katz model. 
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Figure 9. Cell-death effective action cross section including 
effects of nuclear reactions for various ions in C3H10TV2 
cells. 



10-15 [ ... 

100 10 1 102 103 10 4 

Ep, MeV/nucleon 


Figure 10. Cell-transformation effective action cross section 
including effects of nuclear reactions for various ions in 
C3H10TV2 cells. 
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Figure 12. Calculations of proton RBE for survival of 
Chinese hamster cells versus absorbed dose. 


Total dose (with nuclear secondaries) 

Fit to gamma-ray survival curve 
Primary dose 
° 60 Co gamma ray 
a Plateau region- 
1 60-MeV proton 



Figure 14. Survival for V79 Chinese hamster cells irradiated 
by 60 Co gamma rays and plateau-region 1 60-MeV pro- 
tons (from Hall et al. 1978) are compared with calcula- 
tions. 



Figure 13. Depth-dose curve for unmodulated 160-MeV- 
proton beam in water. 


Proton RBE 

I RBE values excluded 
at 95% level of confidence 
(Hall et al. 1978) 

v RBE values unlikely at 

less than 95% level of confidence 
(Hall et al. 1978) 



Figure 15. Proton RBE versus proton dose for V79 Chinese 
hamster cell survival in plateau region of 1 60-MeV- proton 
Bragg curve from first method described in Hall et al. 
(1978). Most likely RBE values fall in space between 
vertical bars. 
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oooo 


Aerobic, theory 

Hypoxic, theory 

Hypoxic, theory, new parameters 

o Aerobic, data (Blakely et al. 1979) 

□ Hypoxic, data (Blakely et al. 1979) 



Dose, cGy Dose, cGy 

(c) E 0 = 395.9 MeV/amu; x = -0.2 cm. (<t) Ei = 395.9 MeV/amu; x = -2.0 cm 

Figure 19. Cell survival as function of dosage at several locations relative to Bragg peak in 2l) Ne beam. 
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